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1. Introduction 

Over the past three years or so the evidence has become overwhelming that 
the cuprate high T c superconductors have a c?-wave pairing state[l]. The 
measurements of Wollman et al.[2] and of Tsuei et al.[3] are especially con- 
vincing since they do not depend on the microscopic physics of the energy 
gap, but instead depend only on the order parameter phase. Other experi- 
ments, such as photoemission[4, 5] and the temperature dependence of the 
penetration depth[6, 7, 8], also strongly support the d-wave picture. 

On the other hand, there is continuing controversy over whether the pairing 
state is a pure d-wave or an s — d mixture [9, 10]. There is indeed evidence for 
a significant s-wave component in YBa2Cu307[ll]. A subdominant s-wave 
component could be compatible with the Wollman et al. and the Tsuei 
et al. experiments provided that it was not too large. The photoemission 
and penetration depth measurements also cannot rule out a small s-wave 
component (either s ± d or s ± id) , although they can possibly put upper 
bounds on the magnitude of the s component. 

Any observations of an s component have important implications for the 
various theories of the pairing mechanism. For example, antiferromagnetic 
spin fluctuations lead to attraction in the d x 2_ y 2 pairing channel, but are 
pair breaking in the s-wave channel[12]. Similarly the Hubbard model with 
a positive on-site interaction U may have a d x 2_ y 2 paired ground state, 
but would presumably not support s-wave Cooper pairs. In the case of 
YBa2Cu307 the orthorhombic crystal symmetry makes some non-zero s- 
wave component inevitable, but a large s-wave component would be dif- 
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ficult to reconcile with either of these pairing mechanisms. On the other 
hand, pairing mechanisms based on electron-phonon interactions, polarons, 
or other non-magnetic excitations (e.g. excitons, acoustic plasmons) could 
be compatible with either s-wave or d-wave pairing states[13]. In these 
models, whether s-wave, d-wave or mixed pairs are more strongly favoured 
would depend on details of the model parameters, and could even vary 
from compound to compound. Indeed there is some evidence that the n- 
type cuprate superconductors are s-wave[14] (or at least they have no zeros 
in the gap |A(k)| on the Fermi surface). This would imply that either the 
pairing mechanism is different for the n- and p-type materials, or that the 
mechanism allows both s-wave or <i-wave ground states depending on the 
band filling. 

In this paper we examine the attractive nearest neighbour Hubbard model, 
which is the simplest model that allows s-wave, d-wave or mixed pairing 
states to occur [15]. We examine the overall phase diagram, paying par- 
ticular attention to the regions near the phase transitions between the s, 
d and s ± id phases. We derive an appropriate Landau-Ginzburg- Wilson 
(LGW) effective action for the model. In the mean-field approximation this 
LGW functional reproduces the usual Hartree-Fock-Gor'kov equations. Be- 
yond mean field theory this functional allows us to examine the large in- 
teraction limit, in which the superconducting phase transition becomes a 
Bose-Einstein condensation of preformed pairs. In this limit there will be a 
pseudogap in the normal state density of states. The strong coupling limit 
is especially interesting in this model because, unlike the case of on-site 
interactions, there are at least two 'species' of preformed pairs, d x 2_ y 2 and 
s. We show below that these become degenerate in the strong coupling 
Bose-Einstein limit. 

2. The Nearest Neighbour Attractive Hubbard Model 

Our starting point is the nearest-neighbour attractive Hubbard model in 
two dimensions, 

H = E ( C U* + H.C) - V^mnj, (1) 

where J2(ij) denotes summation over all of the bonds between the nearest 
neighbour sites. For simplicity we ignore any on-site interaction terms. 

The mean-field gap equations can be used straightforwardly to calculate 
T c for this model (Section 4). The gap equation has stable solutions for 
(extended) s-wave, p-wscve, and ci-wave pairing states, depending on the 
model parameters. One can see in Fig. l(a)-(c) that, for all V, the d-wave 
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Figure 1. Mean-field transition temperature, T c , versus chemical potential, /x, for s-, p- 
and d-wave superconductivity at different interaction strengths; (a) V = It, (b) V = 5t 
and (c) 1/ = 8t. 
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state is strongly favoured at and near to the van Hove singularity at fi = 0, 
while the s-wave state is favoured in the limits of nearly full or empty bands. 
The p-wave solution also has a maximum at \x = but is always small com- 
pared with the d-wave. The mean-field phase diagram also includes a small 
region of s ± id pairing (not shown in Fig. 1) just below the points where 
the d- and s-wave T c cross. These crossing points are therefore tetracritical 
points, where all four phases (normal, s, d and s±id) meet. Below we shall 
closely examine the nature of the phase diagram in the vicinity of these 
tetracritical points. 

If next nearest neighbour hopping terms are also included in the Hamil- 
tonian then there is no longer particle-hole symmetry about half-filling and 
and it is possible to have d-wave pairing predominantly for p-type doping 
and s-wave pairing predominantly for n-type doping. This would be quali- 
tatively consistent with the experimental evidence for s-wave pairing in the 
n-type cuprates[14]. However, for simplicity, below we shall concentrate on 
the case of nearest-neighbour hopping only. 

3. The Landau-Ginzburg- Wilson Effective Action 

In order to examine how fluctuations modify the mean-field picture of Fig. 
1 is is useful to develop a Landau-Ginzburg- Wilson (LGW) effective action. 
This allows us to examine both the weak coupling, V <C 8t, and the strong 
coupling, V 3> 8t, limits on an equal footing. Only in the strong or inter- 
mediate coupling regime will fluctuation effects lead to phenomena such as 
the normal-state pseudogap which is observed in several different high T c 
materials [16]. 

The Landau-Ginzburg- Wilson effective action of the nearest-neighbour Hub- 
bard model can be derived in a way similar to the usual s-wave case of on- 
site interaction[17]. Firstly the grand canonical partition function is written: 



Z 




(2) 
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where Hq is given by the hopping term in the Hamiltonian Eq. 1 (including 
a chemical potential term, —/J.N) and Hj is the interaction term. 



The interacting part may be neatly written in terms of pairing operators 
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on nearest neighbour sites: 



Hi = -y£tr(4.F 4J ), (4) 

(ij) 

where is defined by: 

= ( Cj\ c A c jl c H \ _ (5) 

The physical nature of the interaction is best illustrated by introducing a 
new set of operators, 

(6) 
(7) 
(8) 
(9) 

In terms of these, Fij may be decomposed as 

Fij = ^(Bij-a-Bij), (10) 
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where <r is the vector of Pauli matrices, and 
3 = R 0( 



Bij = Bff, (11) 



Bil s (M^±M, B ,). (12) 

The scalar, £>ij, is even under both time reversal and parity whilst the vec- 
tor, Bjj, is odd under both. We interpret B^ and By as the annihilation 
operators for singlet and triplet Cooper pairs on the bond (ij)[18]. 

In terms of these operators the interaction Hamiltonian becomes: 

Hj = -V£tr(4i^), 

<ij> 

= ~ V H (4^ + 4-%)' ( 13 ) 

<ij> 

where we have used tr(<Tj) = and tr^o^) = 2<5y. 
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Since this is a sum of squared bilinear Fermi operators, the Hubbard- 
Stratonovic transformation [17] can be employed. The Gaussian identity, 

£+ VAU = V_[ ## * e -V-(H 2 +^*+^) ; (M) 
2ni J 

where A is a bilinear Fermi operator and (f> is a c-number, allows us to 
decouple the quartic interaction term in terms of new bosonic fields. In the 
case of Eq. 13 this is accomplished by introducing a complex scalar field, 
ipij(r), and a complex vector field, ^^(r), for each bond in the lattice and 
every imaginary time. The partition function can now be expressed in terms 
of these Bose fields: 

Z = J V{iP,?P*;y,y*}e Sb . (15) 

The effective action is 

S b = -V fdrJ2 (l^(r)| 2 + |^(r)| 2 ) + In J V[c,J]e s f, (16) 
° (ij) 

with the remaining fermions contained in 

S f = f dr[£cl{T)d T c la {T) - H (17) 
J ° 1 To 

~ V H (V£(^(r) + *£(r) • B tJ (r) + H.C) " . 

(ij) 

Equations 15-17 are a formally exact representation of the nearest neigh- 
bour attractive Hubbard model. 

The bosonic fields introduced above are defined separately for each bond 
(ij) in the lattice. It is more convenient to form site-centred combinations 
with a definite symmetry. We define two singlet fields at site r« 

4>d = ^(i>x+^-x-^y-^-y), (18) 

where ip a (^ii t) = ipu+ a (T) for a = ±x, ±y according to the direction of the 
bond (ij). Similarly we can define the symmetrised triplet fields 



(19) 
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In the limit of fields varying slowly in space and time these correspond 
to the Ginzburg-Landau order parameters for (extended) s-wave, d x i_ y i 
and p-wave pairing, respectively. The s- and d- wave pairing fields are even 
under inversion symmetry about lattice site rj, while the triplet fields are 
odd. In the notation of Ref. [19], these fields are the order parameters for 
superconductivity in, respectively, the A± g and B\ g and E u representations 
of the tetragonal point group D^. 



4. The Saddle Point Solutions 

The saddle points of the effective action generate Hartree-Fock-Gor'kov 
mean field theory, where the Bose fields become static and spatially uniform 
order parameters: 



1 

= -pxj E( B «( r *> T )) / > a = x >y ( 21 ) 



where N is the number of sites and (•••)/ denotes self-consistent averaging 
with respect to the fermionic part of the action, Sj. 

The transition temperature for a given order parameter, in the absence 
of any of the others, is given by the solution of 

1 = — V ^ tanh — t- , 22 

in which the weighted density of states, N a (e), is 

^ Q = ^EC Q (k)C a (k)<5(6-6 k ). (23) 
k 

The form factors reflect the point group symmetries of the order parame- 
ters: 

cos^) + cos(ky) a = s, 

C a (k) = t cos ^ ~ cos ( k y) a = d, , 24 , 
sm(k x ) <\ p,.. 

sin(ky) a = p y . 

The solutions of Eq. 22 are shown in Fig. 1 for three values of V. The s- 
wave solution dominates for small and large fillings with d-wave dominant 
near the van Hove peak at the centre of the band. The p-w&ve solutions 
are sub-dominant everywhere. It is expected that interaction with the large 
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d-wave order parameter will suppress the p-wave T c even further. It is clear 
from this that p-wave pairing is irrelevant in the bulk superconductor, and 
henceforth we ignore it and concentrate on s- and d-wave pairs only. 

5. Beyond the Saddle Point: s- and d-wave Mixing 

Starting with the effective action, <S&, we Fourier transform in space and 
(imaginary) time and integrate out the fermions to give a purely bosonic 
action, 

S b = -^ElV^)| 2 + Trm(l-VG ), (25) 

q,a 

where q =(q,iu;), and the trace is over both space-time and spinor indices. 
The fermions of the original theory live on in the form of the Nambu Green's 
function matrix, 

G (k,k>)= ( G ° (/C) _ G ° (jfe) )<W, (26) 

in which Go(k) = (iujn — e^+fi)^ 1 is the Green's function for non-interacting 
fermions. The interaction of fermions and bosons occurs through the po- 
tential matrix, 

w{k ' k ' ] = TP? ( r a (-l + k') M \ k) ) C ^ k " (27) 

where the Einstein summation convention has been used for repeated Greek 
indices and, 

Ck,k' = \ ( cos(A; x ) + cos(*4)) ± \ ( cos(fc w ) + cos(*£)) , a = s,d. (28) 

When ip a is constant in real space (i.e. k = k'), this takes a particularly 
simple form, 

C(k) = Ck,k = cos(^) ± cos(^), a = s, d (29) 
as seen in the weighted densities of states. 

Near T c , where the ip a are small, we expand the logarithm as a power 
series up to fourth order; 

4 , 

Tr ln(l - VG ) = - ^ ~ ^ ( VG o) m + O (V> 5 ) . (30) 
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The odd terms in the series are exactly zero but the even terms survive. 
The quadratic contribution is 

Tr (VG ) 2 = V 2 J2 r a (q)MQ)x a ^(q), (31) 
where the susceptibility is: 

X a(S (q) = -^E G o(fe)GS(fc + (7 )C ik+q Cf +q ,k ) (32) 

' k 

1 \ - ( /(fk+q -») + /(Ck - M) - 1 \ r a r (i 

^ + e k+q + e k -2/x J ^ W" 

The quartic term is 

Tr ( vG °) 4 = ^E^(#(^;(«'')*(9-9'+«v w (^w'). 

(33) 

It is sufficient to evaluate the four body susceptibility at q = q' = q" = 
and then to treat it as a constant, x = x(0> 0, 0): 



^ = ^El G o(fc)| 4 C a (k)C /3 (k)C 7 (k)C 5 (k), (34) 



X 



where £ = e k — /U. 

Thus, after a trivial rescaling of ip by V? , the bosonic action to fourth 
order reads: 

S b « -^E^ + ^^C^^) (35) 



2 

q,af3 



In order to derive a Landau-Ginzburg- Wilson functional from <S&, we expand 
to lowest order in small iui u and | q| : 

2 

(8°* + Vx"P(qj) ~ a^-kf*V + E A". (36) 



10 



where: 



a aP = 

d a ? = 
1 



2<3 



1 _il V /tanh(^/2)^ 

2 ^ 2iV ^ ^ e 



C Q (k)C^(k) I , 



V _ /tanh(^/2) \ _ 



E 



SiV^V £ 2 



V 



167V 



C a (k)c^(k) 



(37) 
(38) 
(39) 



+ 



ld_\ / tanh(/3g/2) \ 
<^ 2 + £#H £ 7 



rc^c^k) 



+2 



+2 



deA d_ Aanh(/3g/2) \ g(Cg, k+q Cf +q , k ) 
9* J < V £ / 

/ tanh(/3g/2) \ d 2 (Cg, k+q C^ +q , k 



g=0 



<9<? 



9=0 



where £ = e k — ^. The last two terms in Eq. 39 are due to the q dependence 
of C a - 



The resulting LGW functional is shown below expressed in real space: 



S LG w = - [ cZr^[/C + V] 
Jo r 



(40) 



where K, is the 'kinetic' part and V the 'potential' part of the action. They 
are given by, 

K = d ss rM q + ^^ + d d U*M d + ^^ (41) 



2m ss ' 2)11,1,1 



V 



where, 



2m sd 

a ss \ip s \ 2 + a dd \i{j d \ 2 + b s \^ + b d W A 



+ k\iPs\ 2 \M 2 + 1 ty.MWd + cc.) 



v 2 

b s = —X ssss , 



V 2 

b d = —X dMd , 



V 2 



ssdd 



(42) 



(43) 



In the static and spatially uniform limit, V may be interpreted as the Lan- 
dau free energy. We will examine the phases arising from this before con- 
sidering further the effects of fluctuations. The free energy given in Eqs. 41 
42 is of the form found earlier by Joynt[20] and by Feder and Kallin[21]. 
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6. Landau Theory 

The phase diagram generated by V is found by simultaneously solving, 

dV dV 



dil> s dip d 



0. (44) 



When k = 0, i.e. there is no coupling between the s- and ci-wave order 
parameters, solving a ss = and a dd = as functions of \i gives the critical 
temperature for each. This is equivalent to solving the linearised gap equa- 
tion, Eq. 22. 



For non-zero k it is possible to have mixed phases. To determine the stable 
phase, two parameters are needed. These are: 



b s 


a dd 


2 




b d 








k{1 + \ 


cos(2A#) 


a dd 


b d 


a ss 



(45) 
(46) 



where A6 is the phase difference between ip s and Vd- Fig- 2 shows the 
stable superconducting phase as a function of these two dimensionless pa- 
rameters. Note that in Fig. 2 the definition of K c is different for k > 0, 
where (A0 = vr/2), and k < 0, where (A0 = 0). 



As Fig. 2 shows, depending on the parameter values s, d, s ± d or s ± id 
phases are possible. The transition from the pure s to the pure d phase 
can occur in a number of ways. For K c > 2, a single first-order transition 
separates the two phases, analogous to a 'spin flop' transition. In this case, 
the 0(4) point (corresponding to B s = l,K c = 2) represents the bicritical 
point where the line of first-order transitions s — > d meets the two normal 
— > superconducting second-order lines. When K c < 2 the transition occurs 
via two second-order phase transitions with an intermediate mixed symme- 
try state, either s±d or s±id. The 0(4) point here represents the meeting 
of the four second-order lines at a tetracritical point. 

Using parameters derived from the nearest neighbour attractive Hubbard 
model, the dashed line in Fig. 2 shows the evolution of (B S ,K C ) for V = 2 
as n changes in the region near the crossing of T* and T^. As \x increases 
we move from the far right hand side of the figure, where T=T* <TJ? and 
d-wave is dominant, to the far left, where T=TJ? <T* and s-wave is dom- 
inant. In these extremes, the dominant <i-wave (s-wave) order parameter 
suppresses the sub-dominant s-wave (d-wave) one even though T<T^(Tf). 
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Eventually, however, a mixed s±id phase arises. This behaviour is seen for 
all V. In particular the mixed s ± d phase is never realised for any values of 
B s and K c derived from the nearest neighbour attractive Hubbard model 
unless an orthorhombic distortion is introduced which is beyond the scope 
of this paper. 

7. The Large V Limit 

In the large V limit it is necessary to retain the thermal and quantum 
fluctuations in the Bose fields derived in section 4. The superconducting 
phase transition then becomes the (Kosterlitz-Thouless in 2-dimensions) 
Bose condensation of these quantum fields. A detailed analysis is beyond 
the scope of the present paper, but it is possible to draw some qualitative 
conclusions based on the effective action derived in sections 3 and 5. 

Fig. 1 shows the evolution of the T c for s, d and p-wave pairing as a function 
of jJL for various values of V/t. The small V limit, represented by Fig. 1(a), 
corresponds to the BCS weak coupling limit in which the pairing is a small 
perturbation on the filled Fermi sea. On the other hand as V is increased 
(Fig 1(b)- (c)) the behaviour changes markedly. One can see that the s-wave 



13 



T c becomes significant even at \i = 0, and that the s and d-wave curves 
become increasingly similar as V is increased. This has a simple physical 
interpretation: in the large V limit the pairs become nearly localised on 
bonds, and the most favourable state is a singlet pair (stable compared to 
triplet by an energy of order t 2 /V). As V increases the effective hopping 
for pairs decreases (also as t 2 /V) and so s- and d-wave pairs become nearly 
degenerate. This degeneracy implies that the Bose condensation of pairs 
in the large V nearest-neighbour Hubbard model is qualitatively different 
from that in the on-site attractive Hubbard model. 

Figures 1(b) and 1(c) also shows that for large V superconductivity occurs 
even when the chemical potential is below the bottom of the electronic band 
at —At. For these values of \i the fermions can only exist as bound pairs, and 
not as free fermions. This is analogous to the large V on-site s-wave case 
discussed by Randeria, Duan and Shieh [22] , in which the criterion for Bose 
condensation to occur in the low density limit was shown to be the same 
as the criterion for formation of a two particle bound state. Figure 1(c) 
implies that in the low density, jjl < —At, limit two-particle bound states 
are formed. An extended s-wave superconducting state then occurs when 
these preformed pairs Bose condense at temperatures below the mean-field 
T c shown in Fig. 1(c). 

8. Conclusions 

We have derived an appropriate Landau-Ginzburg- Wilson effective action 
to describe the nearest neighbour attractive Hubbard model. This allows 
us to discuss both the small V BCS and the large V Bose-Einstein con- 
densation limits on the same footing, as was done for the on-site attractive 
Hubbard model some years ago[23]. The nearest neighbour Hubbard model 
is especially interesting because the phase diagram includes regions of both 
(extended) s-wave superconductivity and ci-wave pairing. We find that p- 
wave pairing is never stable. We have studied closely the region of the phase 
diagram near the cross-over from s- to d-w&ve pairing, and find that the 
two phases are always separated by two second-order phase transitions with 
an intermediate phase of s ± id superconductivity. 

In future it will be interesting to examine more carefully the large V Bose- 
Einstein condensation limit of this model, since the extended s-wave and 
d-w&ve pairing states become nearly degenerate. It will also be interesting 
to see how closely features of this model, such as the existence of s, d and 
s±id phases, correspond to the actual experimental features of the cuprates. 
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